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Abstract: We study nonlinear modes in subwavelength slot waveguides 
created by a nonlinear dielectric slab sandwiched between two metals. We 
present the dispersion diagrams of the families of nonlinear plasmonic 
modes and reveal that the symmetric mode undergoes the symmetry- 
breaking bifurcation with the energy primarily localized near one of the 
interfaces. We also find that the antisymmetric mode may split into two 
brunches giving birth to two families of nonlinear antisymmetric modes. 
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1. Introduction 

Recent progress in nanofabrication opens novel opportunities for engineering smaller optical 
devices. One of the rapidly emerging fields is the study and design of optical integrated circuits, 
which would allow increasing functionalities of basic operating elements in information pro- 
cessing. However there exist some fundamental limits for scaling of optical elements and their 
overcoming is a challenging physical and engineering problem. It is believed that incorporating 



metals, compatible with nowadays electronics, into optical elements would allow overcoming 
fundamental diffraction limits by surface plasmon excitation, squeezing light at subwavelength 
scale, and increasing nonlinear response for switching and signal processing applications. Thus 
plasmonic waveguides and elements are of great interest these days. 

Some basic linear plasmonic elements and devices have been proposed recently JTJ . One 
of the simplest plasmonic waveguides is an interface between metal and insulator supporting 
plasmon polaritons; however, due to losses in metal excited plasmon propagates very short dis- 
tances. Introducing three-layer system helps to increase substantially the propagation distances 
due to the coupling of plasmons at the neighboring interfaces and the field concentration in 
dielectric rather than metal [2 |. Thus, two possible geometries seem interesting for guiding of 
plasmons, namely, insulator-metal-insulator and metal-insulator-metal. In past decades rigor- 
ous linear analysis of these structures has been presented J3]|4]|5|. However, to the best of our 
knowledge, nonlinear waveguides of this type have not been analyzed yet. We can mention only 
a few papers concerning the study of nonlinear plasmons at the metal-dielectric interface where 
the dielectric possesses the Kerr nonlinear response J6] [7] [8] |9) . 

Since only TM electromagnetic waves can be supported by a metal-dielectric interface, the 
corresponding nonlinear Maxwell equations involve two components of the electric field, and 
they can't be solvable analytically in a general case and several approximations and simplifica- 
tions have been employed. The first, uniaxial approximation, proposed by Agranovich et al. J6), 
allows to solve Maxwell's equations when nonlinearity depends on the longitudinal component 
of the electric field. However, in real physical systems the longitudinal component is weaker 
than the transverse component and, therefore, this assumption is valid in specific cases only. 
Stegeman et al. assumed that nonlinearity is caused only by the transverse field component |7); 
this approximation has also some limitations. Numerical studies of the nonlinear problem with 
a single interface has been presented in Refs. J8] |9) while some other studies iflOl [TT1 pre- 
sented the mathematical analysis of the Maxwell equations in quadratures, but the results of 
this analysis are hard to generalize to multilayered systems. 

Based on the uniaxial approximation Q, the studies of long range plasmon-polaritons in 
metallic films embedded into nonlinear materials have been carried out as well lfT2l . But, to the 
best of our knowledge, the modes of nonlinear slot waveguides have been not studied yet. In 
this paper, we provide detailed numerical analysis of the problem of nonlinear plasmonic slot 
waveguide calculating all possible modes in the metal-insulator-metal structure including the 
mode splitting and bifurcations. We demonstrate that the symmetric nonlinear mode undergoes 
bifurcation via symmetry breaking, and a new asymmetric mode with the energy localized near 
one of the interfaces emerges. In addition, we demonstrate an interesting effect of the splitting of 
antisymmetric modes that is associated with the multivalued solutions for linear guided modes, 
which were not discussed in earlier literature. This results suggest interesting applications for 
nonlinear switching of plasmonic guided modes. 

The paper is organized as follows. In Sec. |2]we formulate our problem and discuss the main 
equations and dispersion relations for linear systems. In Sec. [3] we consider the guided wave 
solutions for the corresponding nonlinear problem. Finally, Sec. |4]concludes the paper. 

2. Linear guided modes 

We analyze the guided modes in a symmetric slot waveguide created by a three-layer structure: 
metal - insulator - metal. The schematic of the structure is presented in the inset of Fig. Q] 
For simplicity of the mode analysis, we assume that metal is lossless. This assumption can be 
justified in the case of low attenuation in the frequency range where losses do not change much 
the shape of the transverse modes studied below but affect their propagation length. 

In such structures only TM modes exist, so that we can present the field components in the 
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Fig. 1. Dispersion of the linear guided modes shown as the dependence of the normalized 
guided index vs. normalized slot width, for different values of p = |e„,|/e^: p = 1.3 (top) 
and p = 1.05 (bottom), for symmetric (solid) and two antisymmetric modes (dashed and 
dotted). Inset shows the schematic of the structure, shaded regions correspond to metal. 



form, H = e y H y (x)exp(~jj5z) and E = [e x E x (x) + je z E z (x)]exp(—jl3z)', also we have consid- 
ered that E z is shifted in phase with respect to E x and H y JSJ. Maxwell's equations for TM 
modes can be written in the form, 

dH y n dE- n 

—± = -£E : , pH, = eE x , —±+j5E x = H x , (1) 
dx dx 

where /3 is the effective index (propagation constant), and the coordinates x and z are normalized 
to 2n/X, where A is free space wavelength. 

The variation of the refractive index across the structure is taken in the following form, 

(Ed, xe(0,d) 
\e m ,x£(0,d) 

where d is the thickness of the dielectric layer. 

Implying the boundary conditions of continuity of the field components H y and E z at the 
interface between metals and dielectric, we derive the well-known dispersion relation |4|, 

tanh(X d d) [e„ 2 , Xj + ejA 2 ] + 2e m e d X m X d = (3) 

where X,„ = \/ ft 2 — £,„ and X d = \/ft 2 — 

Detailed analysis of this dispersion relation |3] |4) demonstrate the existence of two types of 
solutions for guided modes, with respect to the structure of the magnetic field, symmetric and 
antisymmetric, also refereed in literature as even and odd modes, respectively. 

A systematic analysis J5) demonstrated that the mode existence and structure are defined 
by the values of the ratio p = le,,,!/^, so that for p < 1 there exists only one (antisymmetric) 
mode with higher cutoff with respect to the slot width d. For p > 1 two modes exist, symmetric, 
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Fig. 2. Dispersion of nonlinear guided modes shown as the total power flow vs. guided 
index for A = 480nm, e,„ = -8.25, e d = 7.84, p = 1.05, a = 1.4 x KT 18 (m 2 /V 2 ), and 
d = 25nm. 

without cutoff, and antisymmetric, with lower cutoff. However, our analysis reveals that at p 
close to 1 and certain slot widths, the lower branch becomes multi-valued so that formally there 
exist three modes for a fixed value of the slot width, one symmetric and two antisymmetric. 

Figure Q] presents the modal dispersion as a function of the slot width d, at different values 
of p. As was predicted earlier Q, at p ^> 1 two modes exist: symmetric (solid curve), and 
antisymmetric (dashed curve). For p approaching 1 we observe three modes, see Fig. Q] (bot- 
tom), symmetric mode and two antisymmetric modes. One of the antisymmetric modes is more 
confined to the interfaces and has lower cutoff (dashed curve), the other mode is less confined 
(dotted curve), and it appears for certain values of the slot width. Since metals are frequency 
dispersive, we can always find a frequency range for different insulators where all those regimes 
can be realized. 

3. Nonlinear guided modes 

We continue our analysis and now assume that the insulator is a nonlinear dielectric with Kerr 
nonlinear response. To be more specific, we consider that the nonlinear dielectric is chalcogen- 
ite glass AsoSe^ with self-focusing nonlinearity sandwiched between two silver slabs. In this 
case, dielectric permittivity of the nonlinear slab can be presented as follows: 

£nln = £,/ + « (4) 

where e d = 7.84, and a = 1.4 x 10~ 18 (m 2 /V 2 ) is the nonlinear coefficient. 

We choose p = 1.05, so that three mode regime is realized, see Fig. Q] (bottom), thus e„, = 
—8.25. For silver this is realized at v = 0.63 x 10 l5 Hz (free space wavelength X = 480nm) lfT4l . 

We solve Eqs. ([TJ numerically by applying the numerical shooting method and find the mode 
profiles for different values of the guided index and nonlinear parameter. To analyze the non- 
linear modes we calculate the total energy flux per unit length in the direction of propagation, 

S = J [[ExH]xz]dx. (5) 

and plot this value for three different values of the slot width, corresponding to different types 
of linear guided waves shown in Fig.Q] 



H, kA/m 



H, kA/m 




Fig. 3. (a-d) Characteristic profiles of nonlinear plasmonic modes shown as the guided 
magnetic field for different branches of the dispersion curves marked by points in Fig. [4] 
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Fig. 4. Dispersion of nonlinear guided modes shown as the total power flow vs. guided 
index for d = 35nm. In the linear limit of vanishing power, the structure supports one 
symmetric and two antisymmetric plasmonic modes. 



First, we analyze the guided modes when there exists only one symmetric mode in the linear 
regime. For this case, we observe that for all values of the guided index /3 > ^/e^ there appears 
an antisymmetric nonlinear mode, see Fig. [2] Even the total flux for this mode may vanish, this 
nonlinear mode has no linear analog, and it requires a finite power to be excited in the structure. 
Our analysis reveals that for smaller slot widths (not shown here) the total energy flux of this 
antisymmetric mode is negative for all guided indices. With the slot width grow, the total energy 
flux of this mode becomes positive for some values of j3, see Fig. [2] The typical mode profile 
is shown in Fig. [2c). 
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Fig. 5. Dispersion of nonlinear guided modes shown as the total power flow vs. guided 
index for d = 50nm. In the linear limit of vanishing power, the structure supports one 
symmetric and one antisymmetric plasmonic modes. 

For j3 ~ 15, a symmetric mode appears, and for larger j3 the flux increases up to its maximum 
(S ~ 6mW /m) at /3 ~ 17 but then decreases again. Also, at j3 ~ 16 we observe a symmetry 
breaking bifurcation (at S ~ 5mW/m) which leads to the appearance of an asymmetric mode 
with power decrease monotonously with /3 . For larger /3 the field becomes strongly confined to 
the interfaces, thus the interaction between plasmons localized at different interfaces becomes 
weaker, so that the dispersion resembles that of a single interface plasmon J9)- Typical mode 
profiles of symmetric and asymmetric modes are presented in Fig. [3ja,b), respectively. 

For larger values of the slot width, three modes appear excited and the antisymmetric mode 
splits into two branches. With larger d (see Fig. [4), the character of the symmetric mode and 
bifurcated asymmetric branch does not change much. However, the bifurcation point is ob- 
served at smaller values of power. Also, the asymmetric mode reaches its maximum. The mode 
profiles are shown in Fig.[3ja-d) for the marked points (1-4). Long-wave antisymmetric mode 
(see Fig.|31 dotted) has a negative flux because the field resides mainly in metal, see Fig. [3d), 
and the mode is less confined to the interfaces in comparison to the short-wave antisymmetric 
mode, [see Fig. [4] dashed; and Fig. Etc)]. We note that for all slot widths there is a range of 
guided indices where several modes at different powers can coexist. 

Further increase of the slot width leads to broadening of the gap between the antisymmet- 
ric branches bringing to the degeneracy of long-wave antisymmetric branch since it becomes 
radiative. The bifurcation power decreases, and it is observed at S ~ ImW/m, see Fig. [5] 

4. Conclusions and acknowledgements 

We have studied the families of plasmonic modes in nonlinear slot waveguides and predicted 
the symmetry-breaking bifurcation of the symmetric modes with the critical power depending 
on the slot width. We have also discussed a complex structure of the asymmetric plasmonic 
modes that original from the splitting of linear plasmonic modes. 
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